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ABSTRACT 

Context. The launch of the Fermi Gamma-ray Space Telescope in 2008 opened new perspectives in the multi-wavelength 
studies of neutron stars, with more than 100 7-ray pulsars now detected. While most Fermi pulsars have been already 
observed in the X-rays with Chandra and XMM-Newton, optical observations with 8m-class telescopes only exist for a 
tiny fraction of them. 

Aims. Here, we aim at searching for optical emission from the Fermi pulsar PSR J1028— 5819 (P — 91.4 ms). With 
a spin-down age r = 92.1 kyr and a rotational energy loss rate E ~ 8.43 x 10 35 erg s _1 , PSR J1028— 5819 can be 
considered a transition object between the young, Vela-like pulsars and the middle-aged ones. At a distance of ~ 2.3 kpc 
and with a relatively low hydrogen column density, PSR J1028— 5819 is a good potential target for 8m-class telescopes. 
Methods. Due to its recent discovery, no optical observations of this pulsar have been reported so far. We used optical 
images taken with the Very Large Telescope ( VLT) in the B and V bands to search for the optical counterpart of PSR 
J1028— 5819 or constrain its optical brightness. At the same time, we used an archival Suzaku observation to confirm 
the preliminary identification of the pulsar's X-ray counterpart obtained by Swift. 

Results. Due to the large uncertainty on the pulsar's radio position and the presence of a bright (V = 13.2) early 
F-type star at & 4" (Star A), we could not detect its counterpart down to flux limits of B ~ 25.4 and V ~ 25.3, the 
deepest obtained so far for PSR J1028— 5819. From the Suzaku observations, we found that the X-ray spectrum of the 
pulsar's candidate counterpart is best-fit by a power-law with spectral index Tx = 1.7 ± 0.2 and an absorption column 
density JVh < 10 21 cm -2 , which would support the proposed X-ray identification. Moreover, we found possible evidence 
for the presence of diffuse emission around the pulsar. If real, and associated with a pulsar wind nebula (PWN), its 
surface brightness and angular extent would be compatible with the expectations for a ~ 100 kyr old pulsar at the PSR 
J1028-5819 distance. 

Conclusions. A far more accurate radio position for PSR J1028— 5819 is necessary to better determine its position 
relative to Star A. Future high-spatial resolution observations with both the HST and Chandra are better suited to 
disentangle the optical emission of PSR J1028— 5819 against the halo of Star A and to confirm the existence of the 
candidate PWN. 

Key words. Optical: stars; neutron stars: individual: PSR J1028— 5819 



1. Introduction 

The launch of NASA's Fermi Gamma-ray Space Telescope 
in June 2008 marked a breakthrough in 7-ray observations 
of pulsars. The Large Area Telescope (LAT; Atwood et al. 
2009) has now detected more than 100 7-ray pulsars, to 
be compared with the 7 detected by the Compton Gamma 
Ray Observatory (CGRO) in 10 years of operation (see, e.g. 
Thompson 2008 for a recent review). Some of them are as- 
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sociated with known, or recently discovered, radio pulsars, 
while others are yet undetected in the radio despite of deep 
searches and are, thus, members of the long-hypothesised 
class of radio-silent 7-ray pulsars, of which Geminga was the 
first example (e.g., Bignami & Caraveo 1996; Abdo et al. 
2010a). In both cases, multi- wavelength follow-ups are im- 
portant to study the pulsar radiation processes over several 
decades in energy, to constrain the source distance and po- 
sition, and to search for pulsar wind nebulae (PWNe), still 
challenging targets in 7-rays despite of the unprecedented 
spatial resolution of the Fermi-LAT (e.g., Ackermann et 
al. 2011). In the X-rays, many Fermi pulsars were already 
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observed prior to their 7-ray detection, while for the oth- 
ers systematic follow-up observations have been carried out 
recently with Chandra, XMM-Newton, and Swift (see, e.g. 
Marelli et al. 2011 for a summary of the X-ray observations 
of Fermi pulsars). In the optical, the source coverage with 
8m-class telescopes is still sparse with only a few Fermi 
pulsars observed beforehand and only a few deep follow-up 
observations on well-selected targets have been carried out 
so far (Mignani et al. 2011; Shearer et al. 2012). 

The radio pulsar PSR J1028-5819 (P = 91.4 ms) was 
discovered by Keith et al. (2008) during a 3.1 GHz sur- 
vey of a sample of low-latitude unidentified EGRET er- 
ror boxes performed with the Parkes radio telescope and 
the Australia Telescope Compact Array (ATCA), prior to 
the launch of Fermi. The pulsar's period derivative P = 
1.61 x 10~ 14 s s _1 yields a spin-down age of 92.1 kyr, a 
spin-down power E ~ 8.43 x 10 35 erg s _1 and a dipolar 
magnetic field B ~ 1.2 x 10 12 G. Thus, PSR J1028-5819 
may be considered a transition object between the young, 
Vela-like pulsars and the middle-aged ones. The dispersion 
measure (DM=96.525 ±0.002 cm -3 pc) implies a distance 
of ^2.3 kpc, according to the Galactic electron density 
model of Cordes & Lazio (2002). Initially associated with 
the EGRET source 3EG 1027-5817 upon its coincidence 
with the 7-ray error box, PSR J1028— 5819 was recently 
identified as a 7-ray pulsar by Fermi (Abdo et al. 2009a) 
and included both in the Fermi catalogue Bright Source 
List (Abdo et al. 2009b) and in the first Fermi-LAT cata- 
logue of 7-ray pulsars (Abdo et al. 2010b). Like for many 
Fermi pulsars, the 7-ray light curve is characterised by a 
double-peak profile, with a phase separation of ~ 0.46, 
which trails the single radio peak by ~ 0.2 in phase. The 
pulsar's 7-ray spectrum is fit by a power-law (PL) with an 
exponential cut-off, with a photon index T 7 = 1.25 and 
a cut-off energy E c = 1.9 GeV. Possible evidence for a 
PWN was also found in the Fermi data (Ackermann et al. 
2011). Recently, PSR J1028-5819 was found to be spa- 
tially coincident with a TeV 7-ray source HESS J1023— 575 
(Abramowski et al. 2011), whose emission is probably as- 
sociated with the PWN. In the X-rays, PSR J1028-5819 
was identified with a faint X-ray source detected by the 
Swift X-ray Telescope (XRT) in a short follow-up observa- 
tion (Abdo et al. 2009a; Marelli et al. 2011) of the 7-ray 
source error box. Unfortunately, the low detection signifi- 
cance (^ 5<t) hampered the characterisation of the pulsar's 
X-ray spectrum. 

While deeper X-ray observations of PSR J1028-5819 
are now in progress with Chandra, no optical observations 
have been performed, so far. Here, we present the first op- 
tical observations of PSR J1028— 5819, performed with the 
Very Large Telescope ( VLT) as a part of a planned sur- 
vey aimed at the optical identification of Fermi pulsars. To 
complement our multi- wavelength analysis, we also present 
a re-analysis of the published SVift/XRT observation of 
the pulsar and the first results of a more recent, yet unpub- 
lished, Suzaku observation available in the public science 
data archive. 

This paper is organised as follows: observations, data 
reduction and analysis are described in Sect. 2, while results 
are presented and discussed in Sect. 3 and 4, respectively. 
Conclusions follow. 



2. Observations and data reduction 

2.1. VLT observa tion 

PSR J1028— 5819 was observed in service mode on 
February 10th and 11th, 2011 with the VLT Antu telescope 
at the ESO Paranal and the FOcal Reducer/low dispersion 
Spectrograph (FORS2), a multi-mode camera for imaging 
and long-slit/multi-object spectroscopy (Appenzeller et al. 
1998). FORS2 was used in imaging mode and equipped 
with its standard MIT detector, a mosaic of two 2kx4k 
CCDs optimised for wavelengths longer than 6000 A. With 
its high-resolution colimator, the detector has a pixel scale 
of 0'.'125 (2x2 binning), which corresponds to a projected 
field-of-view (FOV) of 4.'15 x 4fl5 over the CCD mosaic. 
However, due to vignetting, the effective sky coverage is 
smaller, and is larger for the upper CCD chip. 

Observations were performed with the standard low 
gain, fast read-out mode, with the high-resolution collima- 
tor, and through the high-throughput 6high (A = 4400 A; 
AA = 1035A) and v H igh (A = 5570 A; AA = 1235A) filters. 
To allow for cosmic ray removal, we obtained sequences 
of short exposures (100 s) for a total integration time of 
4800 s in both the &high and which filters. We set such a 
short exposure time to minimise the saturation of a bright 
V = 13.12 star close to the pulsar position, identified in the 
Digitised Sky Survey (DSS) images. Exposures were taken 
in dark time and under photometric conditions, as recorded 
by the ESO ambient condition monitoiQ, with an airmass of 
^1.2 and image quality of ~ 0'.'8, measured directly on the 
images by fitting the full-width at half maximum (FWHM) 
of unsaturated field stars. 

Data were reduced through the ESO FORS2 pipelin^EI 
for bias subtraction and flat-field correction. Per each band, 
we aligned and average-stacked the reduced science images 
using the swarp program (Bertin et al. 2002), applying a 
3(7 filter on the single pixel average to filter out residual hot 
and cold pixels and cosmic ray hits. Since all exposures were 
taken with sub-arcsec image quality, we did not apply any 
selection prior to the image stacking. We applied the pho- 
tometric calibration by using the extinction-corrected night 
zero points computed by the FORS2 pipeline and available 
through the instrument data quality control databasfl 

2.2. Swift observations 

PSR J1028-5819 was pointed by the SVift/XRT on 
November 23rd 2008, during a Target of Opportunity (ToO) 
observation (Obs ID 00031298001) requested soon after the 
pulsar detection by Fermi. The XRT (Burrows et al. 2005) 
is a grazing incidence Wolter I telescope equipped with an 
X-ray CCD imaging spectrometer, yielding a 110 cm 2 effec- 
tive area at 1.5 keV, a 23f6 x 23f6 FOV (2'/36/pixel), and an 
Half Energy Width of 15" in the 0.2-10 keV energy range. 
The integration time was 9.6 ks (7.9 ks accounting for the 
exposure map correction) . The observations was performed 
in PC mode. Data reduction and analysis were performed as 
in Marelli et al. (2011) but using the newest version (6.11) of 
the HEAsof t software package]. On the same visit, the PSR 
J1028— 5819 field was also observed in the near-UV with 

1 http : / / archive . eso . org/ asm/ ambient-server 

2 http : //www . eso . org/ sci/ software/pipelines 

3 www.eso.org/qc 

4 heasarc . gsf c .nasa. gov/docs/ sof tware/lheasof t/ 
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the Ultra- Violet Optical Telescope (UVOT) on board Swift. 
The UVOT (Roming et al. 2005) is a 30 cm, f/12.7 Ritchey- 
Chretien telescope using a microchannel-intensified CCD 
detector which operates in photon counting mode, covering 
the 1600-8000 A range over a 17' x 17' FOV (O'.'5/pixel). 
Observations were performed with the UVM2 (A = 2246 
A; AA = 498A) and the UVW2 (A = 1928 A; AA = 657A) 
filters (Poole et al. 2008) for a total integration time of 4115 
s and 5701 s, respectively. We retrieved the data from the 
Swift science data archive and we combined the single expo- 
sures in order to obtain the best signal-to-noise (S/N). We 
used the UVOT data analysis software available in HEAsof t 
and the version 20110131 of the UVOT calibration files. 

2.3. Suzaku observations 

The field of PSR J1028-5819 was observed also by Suzaku 
(Mitsuda et al. 2007) on July 22nd 2009. The observation 
(Obs ID 504045010) was performed with the pulsar located 
at the nominal aim-point position of the X-ray Imaging 
Spectrometer (XIS; Koyama et al. 2007). The XIS consists 
of three front-illuminated (FI) CCD cameras (XIS0, XIS2 
and XIS3), and one that is back-illuminated (BI; XIS1). 
Each CCD camera is combined with a single X-ray tele- 
scope (Serlemitsos et al. 2007) and has a single chip, yield- 
ing a 330 cm 2 (FI) or 370 cm 2 (BI) effective area at 1.5 
keV, 18' x 18' field of view (l'.'05/pixel), and 1.'8 spatial 
resolution in the 0.2-12 keV energy range. One of the front- 
illuminated CCDs, XIS2, was not available at the time of 
our observation, since it suffered from a damage on 2006 
November 9 and is unusable since then. The observations 
was performed in full-frame mode. 

We downloaded the data from the public Suzaku science 
data archive through the DARTS interfac^E Data reduction 
and analysis of the XIS was performed with HEAsoft fol- 
lowing the procedures described in the Suzaku ABC Guidc0 
version 2. We used the task xispi to convert the unfil- 
tered XIS event files to pulse-invariant channels and then 
xisrepro to obtain cleaned event files. For each XIS, 3x3 
and 5x5 edit modes cleaned event data were combined 
using xis5x5to3x3 and xselect. Following standard prac- 
tices, we excluded times within 436 s of Suzaku passing 
through the South Atlantic Anomaly and we also excluded 
the data when the line of sight was elevated above the 
Earths limb by less than 5°, or is less than 20° from the 
bright-Earth terminator. Moreover, we excluded time win- 
dows during which the spacecraft was passing through the 
low cut-off rigidity of below 6 GV. Finally, we removed 
hot and flickering pixels using sisclean. The nominal on- 
source time is 23 ks. With all the aforementioned data selec- 
tion criteria applied, the resulting total effective integration 
time is 20.2 ks. For the spectral analysis, the response ma- 
trices were generated with xisrmf gen and using ray-tracing 
simulations with xissimarf gen. 

3. Results 

3.1. VLT astrometry 

The radio position of PSR J1028— 5819 has been obtained 
through interferometry observations with the ATCA (Keith 

5 www. darts . isas . jaxa.jp/astro/ suzaku/ 

e heasarc .gsf c .nasa.gov/docs/suzaku/analysis/abc/ 



et al. 2008) and is: a = 10 h 28 m 27!95 (±0?1) and 5 = 
-58°19'05'.'22 (±1'.'5), at MJD=54562. Due to the relatively 
small difference of ~ 3 yrs between the epoch of the refer- 
ence radio position and that of our VLT observation, we 
can neglect any additional uncertainty on the pulsar posi- 
tion due to its unknown proper motion. If assuming the DM 
pulsar distance of 2.3 kpc and a transverse velocity three 
times as large as the peak value of the radio pulsar ve- 
locity distribution (~ 400 km s _1 ; Hobbs et al. 2005), the 
time difference would correspond to an additional uncer- 
tainty (in the radial direction) of ~0'.'3, much smaller than 
the nominal error on the radio position. No more recent 
(and more accurate) radio position for PSR J1028— 5819 
has been obtained by the Fermi Pulsar Timing Consortium 
(Smith et al. 2008), with Weltevrede et al. (2010) assuming 
the original ATCA position of Keith et al. (2008), while the 
measured Swift/XRT position of the pulsar candidate coun- 
terpart, a = 10 h 28 m 27;69 and 5 = -58°19'05'.'16 (Marelli 
et al. 2011), has an estimated radial uncertainty of ~ 5". 
Thus, we assume the ATCA radio position as a reference. 
We note that Simbad reports a quite different position, 
a = 10 h 28 m 39^8 and 5 = -58°17'31'.'0, which is still that 
given in the Fermi-LAT Bright Source List (Abdo et al. 
2009b). 

To register the pulsar position on the FORS2 frames 
as precisely as possible, we re-computed their astromet- 
ric solution which is, by default, based on the coordi- 
nates of the guide star used for the telescope pointing. 
As a reference, we used stars selected from the Guide 
Star Catalogue 2 (GSC-2; Lasker et al. 2008). Although 
the recently released UCAC-3 catalogue (USNO CCD 
Astrograph Catalog; Zacharias et al. 2010) has a better 
astrometric accuracy then the GSC-2, it is cut at brighter 
magnitudes (V Si 16) so that most of the reference UCAC- 
3 stars are saturated even in the short-exposures FORS2 
images. We selected GSC-2 stars evenly distributed in the 
field-of-view but far from the CCD edges and the vignetted 
regions, where geometric distortions appear to be larger 
than expected. We noticed this effect for the first time 
in our data set and, apparently, it affects all FORS2 im- 
ages taken with the high-resolution collimator after 2009 
(M. van den Ancker, priv. comm.). Its origin is still un- 
der investigation. We measured the star centroids through 
Gaussian fitting using the Graphical Astronomy and Image 
Analysis (gaia) tooQ and used the code ASTROn{3 to com- 
pute the pixel-to-sky coordinate transformation through an 
high-order polynomial, which accounts for the CCD distor- 
tions. The rms of the astrometric fits was <r r ~ O'.'l in the 
radial direction. To this value we added in quadrature the 
uncertainty trtr = O'.'l on the registration of the FORS2 im- 
age on the GSC2 reference frames, a tI = \J (3/A s )ctgsC2, 
where ctgsC2 = 0'.'3 is the mean positional error of the GSC2 
coordinates and N s is the number of stars used to compute 
the astrometric solution (Lattanzi et al. 1997). After ac- 
counting for the ~ 0715 accuracy of the link of the GSC2 
coordinates to the International Celestial Reference Frame 
(ICRF), we thus estimate that the overall (la) uncertainty 
of our FORS2 astrometry is 5 T ~ 0'.'2. This is negligible 
compared with the absolute error on the reference pulsar 
position. 



7 star-www. dur . ac . uk/~pdraper/gaia/gaia.html 

8 www . starlink.rl . ac .uk/ star/docs/ sun5 .htx/ sun5 .html 
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Fig. 1. 40" x 40" VLT/FORS2 image of PSR J1028-5819 
in the uhigh filter (4800s). North to the top, East to the left. 
The cross marks the pulsar radio intcrfcrometry position 
(Keith et al. 2008), while the circle marks the position of the 
candidate X-ray counterpart detected by Swift/XKT (Abdo 
et al. 2009a; Marelli et al. 2011). Stars close to the radio 
position are labelled (A-D). The bright Star A (V = 13.12) 
is saturated in the FORS2 image. The intensity scale in 
the image has been stratched to visualise the fainters stars 
(B-D). 



3.2. The search for the pulsar's optical counterpart 

The computed PSR J1028— 5819 position is shown in Fig. 
1, overlaid on the FORS2 V-band image. As mentioned in 
Sectn. 2.1, it falls <> 4" from a bright (V = 13.12) field star 
(Star A) west of it. Unfortunately, since its halo extends 
to 1" from the pulsar position we could not mask Star 
A without risking of severely hampering our observations. 
Indeed, since the angular separation between Star A's halo 
and the pulsar is compatible with the telescope pointing 
accuracy and smaller than the absolute uncertainty on the 
pulsar's position (see above), this could have been acciden- 
tally covered by one of the FORS2 occulting bars (<~ 12"5 
width with the high-resolution collimator). Moreover, Star 
A is part of a triplet (Star B, C) of relatively bright stars 
which was not resolved in the DSS-2 images, the only ones 
available for a quick look of the pulsar field, due to their 
coarse spatial resolution. This would have caused the mask 
to be put even closer to the pulsar position, thus increasing 
the risk of an accidental occultation of our target. We note 
that the northern star of the triplet (Star B) overlaps the 
pulsar position. However, its relatively high flux (B ~ 20.3; 

V ~ 17.5) obviously rules it out as a possible candidate 
counterpart to the pulsar (e.g., Mignani 2011). The same 
argument also applies to the southern star of the triplet 
(Star C) which is of comparable brightness (B ~ 20.7; 

V ~ 17.5). 

A fourth, fainter star (Star D) is visible in the V-band 
image at ~ 2" from the nominal radio pulsar position, su- 
perimposed to the halo of the much brighter Star A, while 



in the B-band image it can not be clearly resolved, par- 
tially because of the slightly worse image quality. The off- 
set is comparable with the la uncertainty on the pulsar 
position, which, in principle, makes it a possible candidate 
counterpart. We computed the flux of Star D through cus- 
tomised aperture photometry using a small aperture of 4 
pixel radius (0"5) to minimise the contribution of the bright 
Star A's halo and optimise the S/N. We sampled the local 
background within circular areas of the same radius as the 
photometry aperture placed below and above the Star D's 
position. We then computed the aperture correction using 
the growth curves of a few unsaturated field stars. The com- 
puted magnitude is V ~ 20.5 ± 0.2, where the photometric 
error is dominated by the uncertainty on the subtraction of 
the background of the bright Star A's halo, which makes it 
also an improbable candidate pulsar counterpart. 

No other object is clearly visible at, or close to, the radio 
pulsar position. Although our images were taken under sub- 
arcsec seeing conditions and with short exposure times to 
minimise the flux contamination of Star A, its bright halo 
still significantly affects the background close to the pulsar 
position, hampering the detection of its faint counterpart. 
We tried to subtract Star A using the model image point 
spread function (PSF), computed by fitting the intensity 
profile of several unsaturated field stars. To this aim, we 
used the tools daof ind, nstar, peak and substar available 
in the IRAF package daophot. Since Star A is partially 
saturated, this procedure obviously leaves some residuals at 
the PSF core, but can account for the subtraction of most 
of the PSF and its wings. Similarly, we also PSF-subtracted 
stars B, C, and D. 

We did not detect any excess signal in the PSF- 
subtracted images, which can be associated with emission 
from a point-like source at the pulsar position. Thus, we 
concluded that PSR J1028— 5819 was not detected in our 
VLT images. We computed the optical flux upper limits in 
the PSF-subtracted B and V bands images from the rms 
of the background sampled in the pulsar radio error box 
within cells of 7 x 7 pixels, i.e. of size equal to the computed 
image FWHM, and we applied the computed aperture cor- 
rection. We then derived 3ct upper limits of V ~ 25.3 and 
B ~ 25.4. Unfortunately, these limits are affected by the 
residuals of the PSF subtraction of the four stars, with Star 
A and B contributing the most, which increase the local sky 
background by a factor of 4 with respect to other regions 
in the field, where the limiting magnitudes are about 1.5 
magnitudes deeper. 

3.3. The pulsar's X-ray counterpart: Swift observation 

It is interesting to compare our derived optical flux up- 
per limits on PSR J1028— 5819, the deepest obtained so 
far, with the extrapolation of the best-fit to the X-ray 
spectrum of the pulsar's counterpart. Unfortunately, this 
is poorly constrained by the available Swift /XRT obser- 
vation (Fig. 2), due to the short integration time (7.9 
ks) and the low count statistics. We extracted all counts 
within a 20"radius using a similar sized background re- 
gion. This yields (2.41 ± 0.69) x 10~ 3 counts s" 1 in the 
0.3-10 keV energy range. We re-analysed the XRT spec- 
tral data using the newest version (12.7) of the xspec 
fitting package (Arnaud et al. 1996). As in Marelli et al. 
(2011), we found a satisfactory fit with a PL with fixed 
photon index Tx = 2, assuming a fixed hydrogen column 
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Fig. 2. Left panel: Swift/XRT image of the field of PSR J1028-5819 in the 0.3-10 keV energy range. The source highlighted 
by the red circle is the pulsar's X-ray counterpart (Abdo et al. 2009a; Marelli et al. 2011). Right panel: Zoom-in on the 
source position. The red cross marks the ATCA position of PSR J1028-5819 (Keith et al. 2008). 



density Ah ~ 5 x 10 21 cm~ 2 inferred from the Galactic 
value along the line of sight using the Heasarc nH too0 
scaled for the pulsar distance of 2.3 kpc assuming a con- 
stant distribution. We note that this value of the Ah is 
about a factor of ~ 3 smaller than assumed by Abdo et al. 
(2009a), Ah = 1.59x 10 22 cm -2 , whose determination, how- 
ever, is not explained. For these spectral parameters, the 
Swift/XRT spectrum corresponds to an unabsorbed 0.5-10 
keV X-ray flux F x = (1.5 ± 0.5) x 10~ 13 erg cur 2 s" 1 . 

3.4. The pulsar's X-ray identification revisited 

The uncertainty on the pulsar's X-ray spectrum obvi- 
ously affects the comparison with our optical flux measure- 
ments. Moreover, the positional coincidence between the 
Swift/XRT error circle of the pulsar candidate counterpart 
and Star A (Fig. 1) suggests that, in principle, this might 
contribute to the observed X-ray emission, which would 
make the characterisation of the pulsar X-ray spectrum 
even more uncertain. As a most extreme possibility, it can- 
not be ruled out that the Swift/XRT source is, actually, the 
X-ray counterpart to Star A and not to the pulsar. The star 
was clearly detected by the Swift/UVOT (see Fig. 3) but, 
admittedly, its possible contribution to the emission of the 
X-ray source was not accounted for by Marelli et al. (2011). 
Thus, we tried to evaluate such contribution by determining 
the star's spectral type via multi-band photometry. 

For the Swift/UVOT observations, we used the FTOOL 
uvotsource to determine the count rate of the star in both 
the UVM2 and in the UVW2 filters. We used a 2"aperture 
in order to minimise contamination from the nearby star 
north-west of it. The count rate was aperture-corrected to 
a radius of 5" in order to apply the count-rate to magnitude 
conversion using the zero points from the the Swift/UVOT 
calibration (Breeveld et al. 2010, 2011). The computed 
magnitudes (in the Vega system) are muvMi — 16.29±0.10 
and mjjvw2 = 16.51 ± 0.11, where the associated errors 



http:/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 



are purely statistic (systematic errors are ~ 0.03 mag- 
nitudes). The star's optical magnitudes are Rp — 12.87 
and V = 13.12, according to the GSC-2 (Lasker et al. 
2008), while the near-infrared magnitudes, as derived from 
2MASS (Skrutskie et al. 2006), are J = 12.38; H = 12.16; 
K s = 12.09. Unfortunately, the UVOT, DSS, and 2MASS, 
images have not enough angular resolution to resolve the 
triplet of stars around Star A, while the star pa 7" north- 
west of it (Fig. 1) is clearly resolved in all of them (see, e.g. 
Fig. 3). Thus, the above magnitudes must be taken with 
due care. However, from our VLT data we see that Star 
B and C are at least a factor of ~ 50 fainter than Star A 
in both the B and V bands. Thus, their contribution only 
marginally affects the estimate of the star colours. 

The colours of Star A are compatible with those of an 
early F-type main sequence star, assuming an interstellar 
extinction corresponding to an E{B— V) ~ 0.09. According 
to the relation of Predhel & Schmitt (1995), this value 
would correspond to an Ah ~ 0.5 x 10 21 cm~ 2 , i.e. lower 
than estimated from scaling the Galactic value along the 
line of sight for the PSR J1028— 5819 distance, suggesting 
that Star A is, possibly, foreground to the pulsar, consis- 
tent with a distance of f=s 1 kpc derived from its distance 
modulus. Usually, F-type stars have an X-ray-to-optical 
flux ratio F x /F v w lO^-lO" 2 (e.g. Krautter et al. 1999) 
which might yield an X-ray flux comparable with that of 
the Swift/XRT source. Thus, it is possible that this is, in- 
deed, the X-ray counterpart to Star A and not to the pulsar. 
However, any conclusion still depends on a better character- 
isation of the X-ray source spectrum and flux. To this aim, 
we have used the archival Suzaku observation (see Sectn. 
2.3). 

3.5. The pulsar's X-ray counterpart: Suzaku observation 

The Suzaku image of the PSR J1028-5819 field obtained 
combining the data from the three XIS cameras is shown 
in Fig. |U Using the sliding-cells source detection algorithm 
in the x image package, a source is detected at the cen- 
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Table 1. Summary of the Suzaku /XIS spectral analysis. Errors are at a la confidence level 
for a single parameter of interest. 



Model 


iV H 
(10 21 cm- 2 ) 


Tx 


kT 
(keV) 


Observed / unabsorbed flux" 
(10" 13 erg cm -2 s" 1 ) 




Power law 
Bremmstrahlung 
Raymond- Smith 

Blackbody 


<1 
<0.5 
<0.6 
<0.5 


1.7 ±0.2 


o + 12 
»-3 

rj+W 

0.5 ±0.1 


1 5 +u - i 

1 -°-0.4 

1 4+ - 1 
0.7 ±0.2 


1.5 ±0.2 
1.4 ±0.2 
1 5+°-i 

1 -°-0.3 

0.7 ±0.1 


1.09/21 
1.16/21 
1.17/21 
1.97/21 



a In the 0.5-10 keV energy range. 




Fig. 3. Same field as in Fig.l but obtained by the 
Swift /WOT through the UVM2 filter (4115 s). The group 
of stars around Star A (see Fig. 1) is unresolved at the 
UVOT spatial resolution. 



tre of the XIS FOV with a signal-to-noise ratio of 18 in 
the 0.2-10 keV energy range. The source best-fit position, 
a = 10 ft 28 m 28f7, S = -58°18'59"3 (J2000), is consistent, 
within the absolute astrometry accuracy of Suzaku (^20"; 
Uchiyama et al. 2008), with the ATCA coordinates of PSR 
J1028-5819 and with those of the source detected by Swift. 

We extracted the source counts from each XIS detec- 
tor using a circular region of 0.'7 radius centred on the best 
fit ximage position. This comparatively small region, corre- 
sponding to ^50% of the encircled energy fraction, was cho- 
sen in order to improve the signal-to-noise ratio of the spec- 
tra and reduce the possible contamination from the diffuse 
emission apparent in Fig. 2) The background was extracted 
from an annulus with radii 3' and 5' centred on the X-ray 
source. Accounting for the counts of all the XIS detectors, 
the 0.5-10 keV net source count rate is (2.7 ± 0.3) x 10~ 3 
counts s , which is not a very high statistics for a de- 
tailed spectral analysis. To increase photon statistics, we 
merged the source spectral files and response matrices us- 
ing addascaspec. The 3-XIS spectrum was binned so as 
that each bin contains a minimum of 15 counts. As done 
for the Swift/XRT spectrum, we performed our spectral 
fits using the XSPEC fitting package (version 12.7). The fit 



with an absorbed PL model (Fig. 5) yields a spectral in- 
dex Tx = 1.7 ± 0.2 and an absorption poorly constrained 
at A H < 10 21 cm -2 at (la), with a reduced xl = 1-09 
for 21 degrees of freedom. This corresponds to an unab- 
sorbed 0.5-10 keV X-ray flux F x = (1.5 ± 0.2) x 10" 13 
erg cm -2 s _1 , very well consistent with the value derived 
by Marelli et al. (2011) from the Swift/XRT observation 
in a similar energy range. We note that our upper limit 
on the N-r is lower than the value inferred by scaling 
the Galactic value along the line of sight (Marelli et al. 
2011), although it is not incompatible with the latter once 
the associated uncertainties are taken into account. We 
also tried other single-component models, like a blackbody 
(BB), a Bremsstrahlung, and a Raymond-Smith plasma 
(Raymond & Smith 1977). However, they all give higher 
X 2 values than the PL, although the fits obtained using 
the Bremsstrahlung and Raymond-Smith plasma models 
are only slightly worse than obtained using the PL (see 
Table Q] for a summary). To improve our results, we tried 
to fit the joint Swift/XRT and Suzaku/XlS spectra, but due 
to the low statistics of the former observation our results 
did not change significantly. 

Our spectral fits confirm the qualitative results of 
Marelli et al. (2011), based on a lower statistics. In par- 
ticular, the good fit with a PL (x 2 , = 1.09) suggests that 
the spectrum of the X-ray source is non-thermal. Thus, it 
would not be produced by coronal emission from a main 
sequence star, like Star A. This would confirm that the X- 
ray source detected by both the Swift/XRT and Suzaku 
is, indeed, the pulsar's X-ray counterpart. Due to the still 
limited statistics, we could not try spectral fits with two- 
component models like, e.g. a PL±BB. Thus, we can not 
rule out that Star A might, indeed, contribute at some level 
to the observed X-ray emission. Deeper observations with 
XMM-Newton will allow one to obtain a more precise char- 
acterisation of the X-ray spectrum of the PSR J1028— 5819 
counterpart, while optical spectroscopy of Star A will pro- 
vide a far more precise classification of its spectral type, 
allowing one to better quantify its expected X-ray flux. In 
addition, high-resolution Chandra images will allow one to 
spatially resolve a possible X-ray source associated with 
Star A from the pulsar's X-ray counterpart and secure its 
identification through the detection of X-ray pulsationj^l 

We also inspected the Suzaku data to search for possible 
evidence of a PWN associated with PSR J1028-5819. As 
mentioned above, diffuse emission on a scale of a few ar- 
cminutes isobserved around the pulsar's X-ray counterpart 



10 We note that the frame times of the Swift/XRT in PC mode 
and of the Suzaku/XIS in full-frame mode are 2.507 s and 8s, 
respectively, which makes it impossible to seach for pulsations 
at the the PSR J1028-5819 91 ms period 
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Fig. 4. Left panel: Suzaku/XIS image of the field of PSR J1028-5819 in the 0.2-10 keV energy range. The pulsar's X-ray 
counterpart is at the centre of the field-of-view. Right panel: Zoom-in on the source position. Like in Fig. 2, the red 
cross marks the ATCA position of PSR J1028-5819 (Keith et al. 2008). 



(Fig. 4), with an apparent elongation in the south-west di- 
rection. The count-rate in this region is (5.3 ± 0.5) x 10~ 4 
counts s _1 arcmin~ 2 in the 0.2-10 keV, computed after 
summing the counts collected with the three XIS detec- 
tors. However, since each XRT/XIS has a different position- 
dependent and broad PSF, with a complicated (roughly 
cross-shaped) and asymmetric profile (see, e.g. Sugizaki et 
al. 2009), it is difficult to determine whether the diffuse 
emission observed around the pulsar's X-ray counterpart is 
real or not. As a test, we compared its surface brightness 
profile with that of a reference X-ray source. Since there are 
no suitable sources detected in the PSR J1028— 5819 field- 
of-view (Fig. 4), we used as a reference an X-ray source de- 
tected in another Suzaku observation (Obs ID 804017020) 
and identified with the ms-pulsar PSR J1231 — 1411. Since 
both sources were observed on-axis, our comparison is not 
affected by the position dependance of the XIS PSF. In 
both cases, we selected photons in the 0.3-10 keV energy 
range in all the three XIS detectors and we used the xissim 
tool to simulate and subtract the background, accounting 
for the detector vignetting. Then, we extracted the radial 
intensity profiles, averaged over a 360° angle, for both the 
PSR J1028— 5819 counterpart and the reference source. 
Finally, we normalised the counts of the reference source, 
computed within a radial distance bin of 10", to that of 
the PSR J1028-5819 counterpart. As seen from Fig. 6, 
the latter seems to feature an excess of counts with respect 
to the reference X-ray source at radial distances ~ 80"- 
200", although such excess is much less evident at larger 
distances. This, coupled with the low count-rate of the dif- 
fuse emission, makes it difficult to determine whether it is 
real, hence possibly associated with a PWN or with an un- 
resolved field source, and to quantify its contribution to the 
observed emission of the PSR J1028— 5819 X-ray counter- 
part, in particular at angular distances smaller than 80". 



10- 



10- ! 
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M+v 







Energy (keV) 

Fig. 5. Suzaku/XIS spectrum of the candidate PSR 
J1028— 5819 X-ray counterpart obtained using data from 
the three XIS CCDs. The background level is marked in 
blue and the best-fit PL model in red. 



Chandra observations will be crucial to confirm the presence 
of the diffuse emission observed by Suzaku. If real, future 
measurements of the PSR J1028— 5819 radio proper mo- 
tion will, then, allow one to search for a possible alignment 
with the semi-major axis of the candidate PWN. 
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Fig. 6. Surface brightness profiles of the PSR J1028-5819 
X-ray counterpart (black) detected in the Suzaku/XIS data 
and of a reference X-ray source (red). The excess with re- 
spect to the reference source is different from zero at the 
~ 5<7 level, (accounting for statistical errors only). This 
suggests the presence of diffuse X-ray emission surrounding 
PSR J1028— 5819 on angular scales of a few arcminutes. 



4. Discussion 

4.1. The pulsar luminosity 

We used our magnitude upper limits in the B-bancQto 
compute the corresponding extinction-corrected optical 
flux and luminosity of PSR J1028— 5819. As a reference 
for the interstellar extinction correction we assumed the 
upper limit on the hydrogen column density (JVh = 10 21 
cm~ 2 ) derived from the fits to the Suzaku X-ray spectrum. 
According to the relation of Predhel & Schmitt (1995), 
this corresponds to an interstellar extinction Ay ~ 0.55. 
Using the extinction coefficients of Fitzpatrick (1999), the 
extinction-corrected B-band upper limit is, thus, Fb Si 
9 x 10~ 16 erg cm~ 2 s _1 . If we assume for the pulsar the 
unabsorbed Suzaku X-ray flux derived from the best-fit 
PL model, F x = (1.5 ± 0.2) x 10~ 13 erg cm" 2 s -1 , our 
extinction-corrected B-band upper limit would correspond 
to an X-ray-to-optical flux ratio Fx/Fb ^ 150, too low to 
independently suggest that the X-ray source detected by 
Swift and Suzaku is indeed the PSR J1028— 5819 counter- 
part. We computed the upper limit on the pulsar optical 
luminosity using as a reference its DM-based distance of 
-2.3 kpc (Keith et al. 2008). We note that while the un- 
certainty on the measured DM is negligible (DM=96.525 
±0.002 cm~ 3 pc), the Galactic electron density model of 
Cordes & Lazio (2002) is known to have an average un- 
certainty of —30% which, however, can be much larger for 
individual objects. As in the first Fermi-LAT catalogue of 
7-ray pulsars (Abdo et al. 2010a), we assumed the typical 
30% uncertainty on the nominal DM-based distance of PSR 
J1028— 5819. This corresponds to an uncertainty of 0.7 
kpc, which we assume in the following discussion. The B- 
band flux upper limit corresponds to an optical luminosity 
Lb £ 9.2 x 10 29 erg s" 1 for PSR J1028-5819, after account- 
ing for the assumed distance uncertainty. Given the pulsar 
spin-down age of 92.1 kyr, we can assume that its, yet unde- 
tected, optical emission be purely non-thermal, as normally 
observed in the youngest pulsars (see, e.g. Mignani 2011). 
Thus, for the pulsar's spin-down power, E — 8.43 x 10 35 

11 For comparison with the luminosities quoted in Zharikov et 
al. (2006), which are computed in the B band 



erg s 1 , our luminosity upper limit implies an optical emis- 
sion efficiency r)B = Lb/E ^ 1.09 x 10 -6 . Both limits are 
consistent with the optical luminosity and efficiency com- 
puted for most pulsars older than —0.1 Myr (Zharikov et 
al. 2006) and confirm an apparent turnover in the pulsar 
optical emission properties. Similarly, the pulsar X-ray flux 
derived from the best-fit PL mode would correspond to an 
X-ray luminosity Lx = (9.5 ± 5.9) x 10 31 erg s" 1 , after ac- 
counting for the X-ray flux and pulsar distance uncertain- 
ties. This value would imply an X-ray emission efficiency 
r]x = Lx/ E = (1.13±0.70) x 10 -4 , for an assumed beaming 
factor / = 1. The derived X-ray efficiency would be consis- 
tent, within the large scatter, with the trend observed for 
the Fermi pulsars (Marelli et al. 2011). 

4.2. The pulsar multi-wavelength spectrum 

We compared the extinction-corrected B and V-band spec- 
tral flux upper limits with the extrapolations in the optical 
domain of the X-ray and 7-ray spectra of PSR J1028— 5819. 
For the X-ray spectrum we used as a reference our best-fit 
PL model with spectral index Fx = 1-7 ± 0.2, while for 
the 7-ray spectrum we assumed a PL with exponential cut- 
off, with spectral index F 7 = 1.25 ± 0.17 and cutoff energy 
E c = 1.9±0.5 GeV, obtained from a re-analysis of the Fermi 
data (see also Marelli et al. 2011). The multi-wavelength 
spectral energy distribution (SED) of PSR J1028-5819 is 
shown in Fig. 7, where we accounted for the lcr uncertainty 
on the extrapolations of the X and 7-ray PL spectra. The 
optical spectral flux upper limits are clearly above the ex- 
trapolation of the 7-ray PL, which does not rule out that 
the optical emission of PSR J1028— 5819 can be associated 
with the low-energy tail of the 7-ray spectrum, as suggested 
for most 7-ray pulsars identified in the optical (Durant et 
al. 2011). This would imply a break in the optical-to-X-ray 
spectrum, as observed in most pulsars (see, e.g. Mignani et 
al. 2010). In this case, the observed optical flux would be 
B = 28-34, a range only marginally reachable by current 
8m-class telescopes. On the other hand, the optical upper 
limits are also compatible with the extrapolation of the X- 
ray PL, which does not allow yet one to confirm the pres- 
ence of a possible spectral break in the optical-to-X-ray 
spectrum. In general, the SED suggests that no single PL 
can fit the optical-to-7-ray spectrum of the pulsar, with at 
least one spectral break required between the X and 7-ray 
energy ranges. This trend has been observed also in other 
Fermi pulsars (e.g., Mignani et al. 2011), probably related 
to a complex particle energy and density distribution in the 
neutron star magnetosphere. 

4.3. The candidate PWN 

Finally, we verified whether the diffuse emission marginally 
detected in the Suzaku data (Sectn. 3.5) would be compati- 
ble, if it were indeed a PWN, with the PSR J1028-5819 en- 
ergetics and distance. PWNe have been detected in the X- 
rays around a number of young and energetic pulsars (see, 
e.g. Kargaltsev & Pavlov 2010 for a recent compilation), 
with luminosities equal to — 10~ 4 -10 -2 of the pulsar spin- 
down power E. In the case of PSR J1028— 5819, one would 
then expect an X-ray luminosity Lx,pwn ~ 8.5 x 10 31 erg 
s" 1 for its PWN, assuming L x .pwn/-E = 10~ 4 . At the PSR 
J1028-5819 distance of - 2.3 ± 0.7 kpc, this would corre- 
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Fig. 7. Upper limits on the extinction-corrected optical 
fluxes of PSR J1028— 5819 compared with the low-energy 
extrapolations of the X-ray (solid black) and 7-ray (solid 
blue) spectral models that best fit the Suzaku and Fermi 
(Abdo et al. 2009a) data. The plotted optical flux upper 
limits are corrected for interstellar extinction based upon 
an JVh = 10 21 cm -2 . The dotted line corresponds to the PL 
X-ray spectrum. The yellow and blue-shaded areas indicate 
the la uncertainty on the extrapolations of the X and 7-ray 
PLs, respectively. 

spond to an integrated flux Fx,pwn ~ (1.3 ± 0.8) x 10 -13 
erg cm~ 2 s _1 for the PWN. By using the WebPIMMS too0 
assuming a PWN spectral index Px,pwn = 2 and an hy- 
drogen column density 7Vh ~ 10 21 cm~ 2 , this would cor- 
respond to a count-rate of (8 ± 5) x 10 -3 counts s _1 in 
the 0.2-10 keV energy range. Following, e.g. Bamba et al. 
(2010), a PWN around a ~ 100 kyr old pulsar would have 
an angular extent of ~ 1U5-21.'5 at a distance of 2.3 ± 0.7 
kpc, i.e. comparable to that of the diffuse emission detected 
in the Suzaku data. For the expected PWN extension, we 
then derive an expected count-rate of ~ (0.2 — 3) x 10~ 5 
counts s _1 arcmin -2 . This value is up to a factor of 30 lower 
than the actual count-rate measured for the diffuse emis- 
sion around PSR J1028— 5819. However, it is not incom- 
patible with it if one assumes, e.g. either a correspondingly 
larger Lx,pwn/E, which would still be within the observed 
range for the known PWNe, or a non-uniform PWN sur- 
face brightness. Thus, it is indeed plausible that the diffuse 
emission detected by Suzaku might be associated with the 
PSR J1028-5819 PWN. 

5. Summary and conclusions 

We performed the first deep observations of the Fermi pul- 
sar PSR J1028-5819 with the VLT. Unfortunately, due to 
the presence of a bright star F-type star (Star A; V = 13.12) 
~ 4" away from the pulsar radio position (Keith et al. 2008) 
we could only set limits of V ~ 25.3 and B ~ 25.4 on its 
optical flux, admittedly affected by the residuals after the 
subtraction of Star A's PSF and by the coarse accuracy of 
the radio coordinates. These limits correspond to an optical 

12 http:/ /heasarc. nasa.gov/Tools/w3pimms. html 



luminosity Lb 9.2 x 10 29 erg s _1 and emission efficiency 
ijb 1.09 x 10~ 6 , consistent with the pulsar spin-down age 
of 92.1 kyr. We also analysed archival Suzaku observations 
of PSR J1028— 5819, which provide a better detection of 
the pulsar's X-ray candidate counterpart and support its 
identification as a pulsar from a better characterisation of 
the X-ray spectrum, which is best-fit by a PL with spectral 
index Tx = 1.7 ±0.2. The inferred pulsar's X-ray efficiency 
(?7X ~ 1.13 x 10 -4 ) would be consistent with the trend 
observed for the Fermi pulsars (Marelli et al. 2011). The 
pulsar multi-wavelength spectrum cannot be fit by single 
PL, implying at least one spectral break in the optical- 
to-7-ray energy range, as observed in other Fermi pulsars. 
We also find a possible evidence for the presence of X-ray 
diffuse emission around PSR J1028— 5819, with a surface 
brightness and angular extent consistent with the observed 
X-ray luminosity and size of known PWNe. More observa- 
tions are required to complete the multi-wavelength study 
of PSR J1028— 5819 at low energies. On the optical side, 
high-spatial resolution observations with the HST, possibly 
in the near-UV to exploit the low interstellar extinction 
along the line of sight, coupled with a much more accu- 
rate radio position, are needed to resolve the pulsar emis- 
sion against the PSF wings of Star A. On the X-ray side, 
follow-up observations with XMM-Newton and Chandra are 
needed to obtain a better characterisation of the pulsar's X- 
ray spectrum, accounting for any possible contribution from 
Star A, and detect X-ray pulsations, which will confirm the 
proposed identification. Chandra observations will also be 
crucial in confirming the existence of the diffuse emission 
marginally detected by Suzaku and study its spectrum and 
morphology. 
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